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Abstract 

The barotropic ideal fluid with step and 5-function discontinuities coupled 
to Einstein's gravity is studied. The discontinuities represent star surfaces and 
thin shells; only non-intersecting discontinuity hypersurfaces are considered. 
No symmetry (like eg. the spherical symmetry) is assumed. The symplectic 
structure as well as the Lagrangian and the Hamiltonian variational princi- 
ples for the system are written down. The dynamics is described completely 



by the fluid variables and the metric on the fixed background manifold. The 

Lagrangian and the Hamiltonian are given in two forms: the volume form, 
which is identical to that corresponding to the smooth system, but employs 
distributions, and the surface form, which is a sum of volume and surface 
integrals and employs only smooth variables. The surface form is completely 
four- or three-covariant (unlike the volume form). The spacelike surfaces of 
time foliations can have a cusp at the surface of discontinuity. Geometrical 
meaning of the surface terms in the Hamiltonian is given. Some of the con- 
straint functions that result from the shell Hamiltonian cannot be smeared 
so as to become differentiable functions on the (unconstrained) phase space. 
Generalization of the formulas to more general fluid is straifgtforward. 
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1 Introduction 



Spherically symmetric thin shells or dust stars (like the Oppenheimer-Snyder one) 
are popular models used extensively in the study of a number of phenomena: prop- 
erties of classical gravitational collapse properties of classical black holes 0], 
quantum gravitational collapse 0, the dynamics of domain walls in early Universe 
@], the back reaction in Hawking effect entropy on black holes or quantum 
theory of black holes 0, |@, to mention just few examples. 

The classical dynamics of objects with discontinuities in matter density is well- 
understood; it is determined by Einstein's equations, the matter dynamical equa- 
tions and some jump conditions at the discontinuity. The jump conditions for the 
step-like discontinuity require that there are coordinates in which the metric is 
at the discontinuity surface 0; for the thin shells, they have been first formulated 
by Dautcourt [T^; Dautcourt's equations have been rewritten in a covariant form 



by Israel [IT 



In many investigations, however, a variation principle, or a Hamiltonian is needed 
from which this classical dynamics follows. Often, such principles (suitable eg. for 
spherically symmetric models) are just guessed from the dynamical equations; some 
attempts to obtain them from more general variational principles are ||12[ and 0. 
Indeed, this is an interesting problem by itself: how the large number of different 
one-degree-of-freedom Hamiltonians scattered in literature is related to the Einstein- 
Hilbert action? For our overal picture of the world has to be self-consistent, even if 
we indulge in using a number of different models, each just applicable for a situation 
under study. 

In the present paper, we reformulate the dynamics of gravitation and ideal discon- 
tinuous fluid in the Hamiltonian form. That is, we identify the canonical variables 
(p's and g's) and Lagrange multipliers and write down a Hamiltonian functional of 
these variables; we show that the constraints and the canonical equations resulting 
from this Hamiltonian are equivalent to the system of Einstein equations and the 
ideal fluid dynamical equations (plus the Israel equations in the case of a thin shell). 

To identify the suitable symplectic structure and find the variational formulas, we 
employ the methods described in detail in |jl3| and their application to general rela- 
tivity as given in |T^. We will, however, keep the paper self-contained by motivating 
and explicitly performing all relevant derivations. 

The model of matter used extensively in this paper is that of the simplest kind: 
the barotropic ideal fluid. This can be formulated as a Lagrangian field theory 
without any constraints |T5[. Generalization to ideal fluid with internal degrees of 
freedom (such as |]I6[) or to any conservative continuum should be straightforward; 
in any case, the gravitational parts of our Lagrangians and Hamiltonians (which 
represent the solution to the main problem) have general validity. 
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In each particular case, the classical dynamics can be obtained from a variational 
principle that has the same form as the corresponding variational principle for a 
smooth system, if some particular generalized functions are allowed to describe the 
matter distribution: the step function for star boundaries and the (5-function for thin 
shells (cf. This simplicity is, however, traded for the freedom in the choice of 

coordinates: the generalized function approach works only if the metric is for 
the step, and C° for the 5-discontinuity. We transform, therefore, the Lagrangians 
and the Hamiltonians to the so-called surface form containing only smooth variables; 
such Lagrangians and Hamiltonians as well as symplectic forms decompose into sums 
of volume and surface integrals. The transformation can best be done in the so-called 
adapted coordinates; these are coordinates in which the embedding functions of the 
surfaces of discontinuity acquire the simplest possible form. The result, however, 
is covariant in the sense that arbitrary smooth coordinates can be chosen inside of 
each separated volume (left or right to the discontinuity surface) as well as along 
the discontinuity surface itself. 

An important trick is used throughout the paper: we work in coordinate systems 
which are always adapted to the position of the discontinuity surface. This way 
the discontinuity surface may be considered as a fixed submanifold of the space- 
time. Thus, the dynamics of the star surface or thin shell is not described by the 
spacetime coordinates of these objects but by the evolution of the physical fields 
like metric of matter fields along the surfaces. Then, for example the variations 
and time derivatives of the embedding functions of the two-surfaces of discontinu- 
ity in the three-surfaces of constant time both vanish identically. Our formulas are 
written only for one hypersurface of discontinuity; an extension to arbitrary many 
hypersurfaces is easy if they do not intersect each other. 

Two interesting problems arise. First, we do not show that the dynamics makes 
sense even on-shell. By that, we mean that there is to be a well-posed initial value 
problem. One ought to be able to define some nice space of initial data, consist- 
ing of those values of the canonical variables that satisfy some well-defined set of 
constraints, jump and fall-off conditions so that a unique solution to the dynamical 
equations will exist in a neighbourhood of the initial surface. In this paper, we 
shall just assume that the dynamics is all right. At least in some special cases (like 
spherical symmetry) , the space of classical solutions is well-known and it is as large 
as one expects. 

The second problem is to show that the Hamiltonian formalism defines a (reg- 
ular) constrained system. This means that one can find a phase space (possibly 
an extension of ours), a complete set of constraints, and a Hamiltonian satisfying 
the following conditions: 1) the constraints and Hamiltonian must be differentiable 
functions on the phase space so that their Poisson brackets are well-defined and 2) 
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the Hamiltonian must be first class and the constraint set must be spht nicely into 
the first and second class constraints (Bergmann-Dirac analysis, cf. Of course, 

such an 'off-shell' formulation is necessary as a starting point for Dirac quantization. 
The difficulty is that some constraints at the shell are not different iable functions 
on the phase space even if they are smeared along the shell, because the smearing 
is then only two-dimensional, whereas the differentiability would require a three- 
dimensional smearing. Without an off-sell formulation, the way to quantum theory 
need not be barred however. One can try to solve the singular constraints and to 
substitute the solution back into the action so that a variational principle results 
which leads to equivalent dynamics without the singular constraints [|1^]. A problem 
with such a procedure seems to be that equations quickly grow very messy. 

The plan of the paper is as follows. Sec. 2 is devoted to the step. Sec. 3 to 
the 5-function discontinuities. Sec. 2.1 introduces the ideal fluid model and its 
dynamics in a fixed spacetime (metric) background. Basic formulas of the Lagrange 
and Euler pictures concerning Lagrangians, Hamiltonians, strees-energy tensors and 
equations of motion are derived; these equations apply to both step and 6 function 
discontinuity. The method of variation formulas is presented, which enables us to 
find the symplectic structure as well as to generate the equations of motion. The 
surface of discontinuity can be moved without problems as far as the metric is 
fixed. In Sec. 2.2, the fluid is coupled to the dynamical gravity. Relevant formulas 
concerning the variation of the Einstein-Hilbert action are collected. The surfaces 
of discontinuity are now fixed. This helps to avoid some formal problems. The 
variation formulas for the system are written in Lagrangian and Hamiltonian form. 

In Sec. 3.1, an action for the thin shell and dynamical gravity is written down 
in the Lagrangian formalism; the shells are fixed and generalized functions are em- 
ployed. In Sec. 3.2, the adapted coordinates are used to transform the (Lagrange 
formalism) action into a sum of volume and surface integrals disposing of the gen- 
eralized functions and gaining more coordinate (gauge) freedom (four-covariance): 
arbitrary coordinates can be chosen left to the shell, right to the shell, and along the 
shell. In Sec. 3.3, the variation of the action in the surface form is calculated and 
the obtained dynamical equations are listed; they contain Israel's equation. The 
variation formula is derived; this is only three-covariant: the foliation by spacelike 
surfaces t = const must be such that the t-surfaces are continuous but can have a 
cusp at the discontinuity surface; the embedding functions of the discontinuity two- 
surface in the t-surfaces must be time- independent. Sec. 3.4 contains a Legendre 
transformation to a Hamiltonian formalism; the general form of the Hamiltonian for 
the system of thin shells and gravity is presented. In sec. 3.5, the explicit functional 
dependence of the Hamiltonian on the dynamical variables is written down and the 
geometrical meaning of the surface terms in the Hamiltonian is disclosed. In Sec. 
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3.6, the variation of the Hamiltonian is exphcitly calculated so that all canonical 
equations and constraints following from the Hamiltonian can be listed. This not 
only enables us to check that the Hamiltonian generates the desired dynamics (in- 
cluding Israel's equation) but also to classify the resulting equations into 'canonical 
equations' and 'constraints.' For example, the six relations that are equivalent to 
Israel's equation consist of one super-Hamiltonian constraint, two supermomentum 
constraints, two singular constraints (these cannot be made differentiable by smear- 
ing), and one canonical equation. Some, necessarily preliminary, discussion of the 
result is given. 

2 Fluid with a step discontinuity 

Our point of departure in this section is the description of relativistic barotropic 
perfect fluid as given in Ref. ||T5[ (observe that this description is easily extended 
to any conservative continuum). We will extend and modify the method so that it 
allows for discontinuous matter distributions admitting such situation like a jump 
of density at the boundary of a star (a step-function type of discontinuity along a 
timelike hypersurface). 

2.1 Fluid in gravitational field 
2.1.1 The description of the fluid 

The fluid that have just 'mechanical' degrees of freedom consists of identifiable 
elementary volumes — mass points of the fluid. It can, therefore, be completely 
described by specifying the mass and the spacetime coordinates of each of these 
mass points. All mass points form the so-called matter space Z, which is a three 
dimensional manifold; let 2;", a = 1, 2, 3, be some coordinates in Z . Let us denote 
the spacetime by M and let x^, = 0, 1, 2, 3 be some coordinates in M. The state 
of the fluid can then be described by a map : M i-^ Z, in coordinates z°-{x^\ the 
particle trajectories are then determined by z^-i^x) = const. The matter space Z is 
equipped with a scalar density h{z), which determines the mole or particle density 
of the fluid, so that the number N of particles or moles in the volume (Z Z is given 
by N(V^) = Jy^ d^z h. We assume further that h has a step discontinuity at a two 
surface in Z, defined by the equation F{z) = 0, where F is a smooth function 
with non-zero gradient Fa- Let S := (^"^(S^) be a timelike three-surface separating 
M in two open subsets and V~ so that h{z{x)) > for x G V~ and h{z{x)) = 
for X G V~^. One can make much more general assumptions (e.g. allowing for several 
matter filled regions) , but this will only complicate the description without requiring 
any new method of approach. 
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The map ( and the density h define mole (particle) current j'^ in M by 



where 

dz"" 



j^{x) is discontinuous at S, j^^ 7^ in V~ , = in V'^ and j^i is tangential to S. 
(We denote the limits to S from inside by the index — .) j'^ is a vector density; it is 
easy to show that j'^ is identically conserved everywhere in M, j^^ = 0. 

The current j^(a;) defines the spacetime four- velocity u^{x) and the rest mole 
(particle) scalar density n of the fiuid in V~ and at S by 

f = ^\nu^, (2) 

where g := det{g^^) and g^yU^u'' = — 1. Hence, 

1 



n = ^J-g^.uj'^f] (3) 

V 1^1 

n has a discontinuity of a step type at S. 

In [|^, it is shown that the fiuid equations of motion can be obtained from the 



Lagrange density which is given by 

Lra = -^/\g\ne{n), (4) 

where e(n) is the energy per mole in the rest frame of the fiuid and is considered 
as a function of z", z'^ and (7^1,. As the specific volume V (i.e. volume of one mole 
in the rest frame) is 1/n, we obtain for the presure p of the fiuid 

p=-— = ne (5) 

in V-. 

2.1.2 Stress energy density 

By definition, the stress-energy tensor density of ideal fiuid (see, eg. 0) has the 
form 

T>^'' = s/\^\{{p+p)u^u-'+pgn, (6) 

where 

p = ne{n) (7) 

is the rest mass density; T^'^ has a step discontinuity at S. In this section, we collect 
some important formulas valid for this tensor density. 
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Let us vary the action of the fluid 

d'^x \ \g\ne{n) (8) 
V- ^ 

with respect to g^^. Using Eq. (|^), we have 

^(V^^) = ^J^gi^^ = )=n-'rfSg^,, 

which, together with the Eq. (|^) and the well-known variation formula for determi- 
nats yields 

6n = -^n{g^'' + u^u'')Sg^.. (9) 

Then, 

61^ = d^x^\ (-^g^'ne + ]^n{g^^ + u^^u^){e + ne')) 
and a straightforward calculation using Eqs. (0), (|^) and (P) leads to 

T^-[x) = 2 , ■ (10) 

The next important relation is the Belinfante-Rosenfeld theorem (|ll9||, pO| and 



2T|| ) applied to our case: the Lagrange density must satisfy the following identity 



This equation is equivalent to the requirement that Lm is a scalar density, and its 
derivation is straightforward. From the identity (p!!]) and the formula (|10|), we obtain 
immediately that 

Fit 

Tti = L^Si: - ^zt. (12) 

Thus, the so called canonical stress-energy tensor density on the right hand side is 
equal to the source of gravitational field. 

The formulas (0) and (|12D imply the Noether identity: 

There are counterparts to Eqs. (|T0D-(|T3D within any description of any type of ideal 
fluid. Derivation of Eq. (|13|) starts from the equation 



If one substitutes for Tj^ from Eq. (|T2|) into the first term on the right hand side 
and from Eq. ( |T0| ) into the second one, the identity follows. One consequence of 
the Noether identity is that the four components of the covector V^T^^ are not 
independent: 

fV.Ti: = 0, 

because the definition of implies the identity j^z'^ = 0. Hence, the equation 
system V ^T^^ = contains only three independent equations (Euler equation); the 
energy conservation equation, 

for the fluid is satisfied identically within our description. 
2.1.3 The variational formula 

Let us consider the 4-dimensional volume V enclosed between two Cauchy surfaces 
and the boundary S of the fluid divides V into two parts, V~ and V~^, and 
similarly Si into S~ and 5*^^, i = 1,2. We assume that Si and 5*2 are C^-surfaces, 
that is the induced metric on, as well as the unit normal vector to. Si and 5*2 are 
both C\ 

Let us vary the matter action with respect to z°'{x) and g^j,^; we obtain 

The first two integrals can be transformed to surface integrals along S, 5*1 and 
5*2. For this aim, we use the coordinates in M that are adapted to the surfaces. 
This means that = ti along 5*^, i = 1,2 and 

= x^ (15) 

along E, so that x'', k = 1,2,3, are coordinates on Si and x", a = 0,1,2 are 
coordinates on S. Then, the change of F if we vary z"" is given by 

F{z%x) + Sz^'ix)) =x^ + FaSz'', 

and the coordinate x^ of S changes by 

Sx^ = -Fa6zX^=o. 
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Thus, we obtain for the first term: 



I5V- 

For the second term, we have 

/ d'xdJ^dzA = I dJ:^6z^+ I dS^6z'^- I dS^5z\ 

Jv- \ OZ'^ J Jt, OZ^ JS2 OZq J Si OZq 

where the abbreviations rfS = dx^dx^dx^ and dS = dx^dx'^dx^ are used. Eq. (|15D 
imphes that FaZ'^ = 6^, so we can write 

7-1 h 

Collecting these results, we obtain the equation 



Thus, the field equations consist of volume equations that hold in V : 



and surface equations that hold at S: 

fdL 



■'m b 



z: - L^6' = 0. 



The surface of the star is an observer independent dynamical element of the system. 

Let us discuss the meaning of the field equations. For the volume equation (|1^, 
we just invoke Noether's identity, Eq. (p!3|); we can then see that they are equivalent 
to the conservation equations V^^T^ = 0. The surface equations can be rewritten as 
follows. First, using Eq. ([I^), we have 



Hence 



Tl^Fbzl = + L^S'}j Fbzl (19) 
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However, F^z^^ is covector normal to S, so the three surface equations (|18|) can be 



written in a covariant form as 

^^1^ = 0, (20) 



where = T'j^rhv and rfiy is any normal covector to S. Of these four equations, 
only three are independent, because T^j'^ is identically zero, as one easily verifies; 
in fact, the three equations T^|s = imply (PU]). It follows further from Eq. 
and from F^iZ^^u^ = that Eq. (|20|) is equivalent to the condition that the pressure 
vanishes, p = 0, at the surface. 

Formula (0) is only valid in the adapted coordinates; in particular, the La- 
grangian density Lm must be expressed in these coordinates. Let us pass to more 
general coordinates. In fact, the volume integrals are already in a covariant form, so 
we have just to transform the boundary integrals. However, we will need explicitly 
only the integrals over the Cauchy surfaces. 

Let us define 

--If- 

The reader can easily verify that Eq. (|21|), which is written in the adapted coordi- 
nates, defines a three-density Pa along Si (independent of the adapted coordinates), 
because the quantity Lm is a four-density. 

The field equations (both volume and surface) are, therefore, equivalent to the 
formula 
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f d'xLm= f d^ypjz'^-f d'ypjz'^ + lf d^xT'^'^Sg^,, 
Jv- J 3-2 J Si 2 Jv- 



and this formula is valid in any coordinates x^ in V~ and y^, k = 1, 2, 3, along Si. 
Let us denote the matter occupied part of the Cauchy surface x° = t by St and 
the corresponding part of the matter space by Z" . The following coordinates will 
simplify all calculations: the intersection fl S is given by x^ = x^{t) and St by 
^3 ^ x^{t); x^ , k = 1,2, 3, are coordinates on Sf, we will call them 'time dependent 
adapted coordinates.' Then, we can write 



d\ Lm= dt f ^ / d^xpaSz^ + 1 I d^x T^'dg^, 

V- Jti \dtJx^<x^{t) 2Jx^<x^{t) 

If we define the Lagragian Cm by 

Cm / dj X 

JSt 

and go to the limit S2 Si, we obtain the variation formula for the Lagrangian: 

6Cm= f d'xipadz'^y + [ d^xpjz'^x' + l [ d^xT^'^Sg^,. (22) 
J St Jast 2 J St 



■'m 
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The variation formula has been derived by careful inclusion of all 'boundary terms'; 
this will be the main strategy for our derivation of the Lagrangian and Hamiltonian 
formalism that admit discontinuities. The role of the variational formula is to gen- 
erate all dynamical equations (including the definition of momenta): the variation 
of the L. H. S. is to be calculated and compared with the R. H. S. 

Eq. (p2D is also the point of departure for the transformation to the Hamiltonian 
formalism. 



2.1.4 The Hamiltonian formalism 



Let us disregard the surface term in the formula (]22|) and define the Hamiltonian 
density Hm by a Legendre transformation of the form 

Then we obtain for the variation of the Hamiltonian Timi which is defined by 

JSt 

the relation 

6Hm = 5 d^xpaZ^-d d^xLm- 
JSt Jst 

Performing carefully the variation in the first term and substituting from Eq. 
for the second, we have 



d^x {dpaz" - PaSz") + [ d^xpaiz^dx^ - x^dz") 
-If d'xT'^''6g,,. (23) 

Z J St 

From this equation, not only the canonical equations are obtained, which will con- 
sist of volume and of surface equations, but also the symplectic structure of the 
system can be read off, which is given by the fist two integrals on the R.H.S.; these 
integals can be interpreted as uj{6z"', 6pa', z"',Pa), where u is the symplectic form and 
(5z'^(x), 6pa{x)), {z'' (x) , Pa{x)) are two vectors; notice that u has a surface part. This 
is a general and very important observation, which will, for instance, help to decide 
what are the canonical variables for the shell in the second part of the paper. The 
last integral in (^) represents variation with respect to non-dynamical 'parameters' 

We will need the following relations between the Hamiltonian density and the 
stress-energy density: 

T° = -H^, (24) 

T'k = -Pazl (25) 




Tt = r-i^Pa-H^\5^+^-^zl (26) 
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which are vahd in the time dependent adapted coordinates. The first two equations 
are obtained immediately from Eq. (0) and the definition of Pa- To derive the last 
equation, we first notice that Eq. (P5|) has the following consequence 



T^^ = -2^. (27) 



Second, we derive an equation analogous to (|TT|) for Hm', we use the fact that Hm 
behaves as a three-density if we change the coordinates keeping t fixed and that 
Hm is a function of z", z^, pa and g^^, = Hmiz"", z^,pa, g^.^)- 

dzf dpa dgpi 

Then, Eq. (^) follows immediately. Again, analogous of Eqs. (p4D-(p7D are valid 
for all types of ideal fiuid. 

To finish the Legendre transformation, we have also to express the velocity in 
terms of Pa, z'^ and z^. in the Hamiltonian. This is not completely straightforward. 
To begin with, we substitute for Lm from Eq. (^ into the definition ( pi]) of pa- 

I dn 



where 



Eqs. (0) and (|^) imply 



dn 



dn 1 dj'^ 

--U, 



dz" f\^\ ^dz"' 



The definition of implies the following identity (cf. [jl5| 



^4=f5:-f5i:. (28) 



Combining the three equations, we easily find 



= (29) 

pr 



where j° depends only on and 2;^ (cf. (|l])): 

f = h{z)det{z'^). 
The (3+l)-decomposition of the metric (see e.g. ^ 

gu = Qkh g = -qN 



k, 

2 
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gives, with the help of Eq. 




) 



2 



q^'ukUi = -1 - 



Substituting for Uk from Eq. (0), we obtain the identity 




) 



2 




1 + 



(30) 



This equation determines n as a function of pa, z"' and 2^. The solution depends on 
the unknown function p{n) and is determined only implicitly, in general. 
The identity (^8]) implies that 



where is the matrix inverse to z'^. As depends linearly on i", we have 



According to the formula (0), this determines the form of the Hamiltonian Tim- 
2.1.5 The Euler picture 

At this stage, we have derived all important formulas of the Hamiltonian formalism 
in which the fluid is described by the functions z"'{x)] these are Lagrange coordinates 
(cf. |]16|) and we can call the formalism 'Lagrange picture\ Sometimes the Euler 



= -3'x\z\ 



or 





(31) 
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picture is more practical, however. This can be obtained by the following canonical 
transformation. The new fields x''{z,t) (Euler coordinates) are defined by 

x''{z{x,t),t) = x'', Va;^^ (32) 

and the conjugate momenta Pk by 

Pk{z, t) := -X{z, t)zt{x{z, t),t)pa{x{z, t),t), (33) 

where 

One easily checks that Eqs. (|32D and (^) define a canonical transformation. 



Let us first derive some useful relations. By differentiating Eq. (p2D with respect 
to t at constant x^, we obtain 

= (34) 

where 

u dx^ 
x„ — 



dz°- 

The derivative of the same equation with respect to x^ at constant t gives 

<4 = Sf. (35) 

If the field z'^{x, t) is changed to 5"(x, t), then x''{z, t) is changed to £^'(2;, t) satisfying 

x''{z{x,t),t) = x^ , Vx^,t. 

Thus, if 5°(x, t) = z°-{x,^) + fc'^(x, t), the above equation implies that 

b^x^ = -x'^Jz". (36) 

The symbol 5* is to stress and to remind us that this variation is of different kind 
than S, if applied to fields: the former is obtained by comparing the values of the field 
at the same point of the matter space, that is at different points of the spacetime; 
the latter compares the values of the field at the same point of the spacetime. 

With the help of the above relations, we can transform all formulas of the Hamil- 
tonian formalism. Let us start with Eq. (|23|) . First, the inverse transformation for 
the momenta follows from Eq. (0): 

Pa{x, t) = -X-\x, t)x^(2;(x, t),t)Pk{z{x, t)t), (37) 

The time derivative of this equation at constant x^ can be calculated with the result 

Pa = [x'^xTzLPkx' + z\xl,Pkx' + x'iZ^Pud.x' + X^zf^bPfc x' 

- Pidax' - xiPu) , (38) 

13 



where we introduced the abbreviation 

and similarly x^^. Analogous formula holds for the variation 5pai one just have 
to replace dots by (5's. Employing these equations, we obtain after a lenghty but 
straightforward calculation 

5paZ'' - PaSz^ = X-'iS^Pkx' - PkSy) + X-'da{PkZtSyx') - X'^ daiPlztS^x" x') , 

or 

/ d^xiSpaZ'^-paSz'') = [ dhid^Pki'-PJ^x") 
JSt JZ- 

+ \ _d''z da ({Pkzt - Pizt)5,xH') . (39) 
Then, we transform the second integral on the R.H.S. of Eq. (p^): 

d'^xpaiz'^dx^ -x^6z^) = I d^xdi(pa{z''5,x^ -x^dz")) 

x3=x3{t) JSt ^ ' 



d^zXz^db (X'^Pk{x^5,x^ - x^5. 



we have used Eqs. (|3^),(p^),(p5|) and ([371). Because of the identity db{Xz\) = 0, we 
obtain finally 

/ d^xpaiz^Sx^ - x^6z^) = - f d^z db UPkz\ - Pizl)5,x^x^) . (40) 

Jx^=x^{t) J ^ ^ 

Eqs. (0), (0) and (EO) imply 



dUm = dh {6,Pk x^ - P^Kx^) + lJ^^ dh XT^^dg^^ixix, t),t). (41) 

Thus, the symplectic form has no surface term in the Euler picture. The variation 
of the metric in the last term on the R.H.S. is independent of the other variations, 
and it is defined by comparing values of the metric at the same spacetime points. 

Let us suppose that Sgfj,i,{x) = 0, and let us introduce the transformed Hamilto- 
nian density H'^ by H'^ = XHm, so that 

f d^zH'^. 

•J z 

Then, 

5Hrr,= [ dh6,H'^, 

Jz- 
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and we have 



Jz- Jz- 



(42) 



In this form, the variational formula is suitable for derivation of the canonical equa- 
tions. To this aim, let us calculate 5^H'^] H'^ is of the form H!^{x^ , x^, Pk) , hence 



'J a' 



dx^ 



dx^ 



f)H' 

S.X" + -T^6,Pk + da 



'dH' 



dx'l 



Thus, the field equations consist of the volume equations: 



x' 



-Pk 



and the surface equations: 



dPk' 

dx^ dx^ 
dH' 



dx^ 



0. 



Let us check that Eq. (^Sf ) is equivalent to (pO|). We have 

H'^ = XHUz^ix'), zt{xi),p,{z%x'), 4(4), Pk),g,Ax')), 

so that 



dx^ 



\P — ^^kHm + ^ 



dz\ '^'^(9x^ ~'~ dp, 



dxl' 



Eqs. (^) and (|37|) imply 



d± 

dxl 
dpb 
dx^ 



and we obtain easily 



d_Kn 

dx^ 



Xz? H„ 



dH^ \ H dH^ A 



dH, 



dp, 



dz\ 



dpa 



m b 
PbZk- 



Application of Eqs. (|25| ) and ( ^61) as well as (|23D simplify the expression to 



dH^n 
dx^. 



-^z^Ij^ . 



(43) 
(44) 

(45) 



(46) 



Hence, Eq. (^51) becomes 
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which is equivalent to (|2y) . 

Finally, the transformed Eq. ( |30D reads 



For example, in the case of dust, p = pn, where /i is a constant (rest mass per mole 
or particle) and Eq. ( ^Tf ) can be solved explicitly: 

1 fih^ 

n 



For dust, Eq. (^I]), which determines the form of the matter Hamiltonian, specializes 
to 



T° = -N^Jp^ify + q^hlz\paPb - N'ztpa. (48) 
2.2 The gravity becomes dynamical 

In the foregoing sections, gravity was just an external field. Here, it will become 
dynamical: the metric g^u{x) will satisfy Einstein's equations with the fluid stress 
energy tensor as a source. 

2.2.1 Description of the system 

The main problem which we shall meet is the following. If Einstein's equations are 
satisfied, the discontinuity in the distribution of the fluid leads to a discontinuity in 
derivatives of the metric. Thus, we must allow for such discontinuity from the very 
beginning. Moreover, a general variation of the metric, which includes a shift of the 
coordinates of the discontinuity, will have a jump of higher order than the metric 
itself: if the second derivatives of the metric have a jump, then the first derivative 
of its variation will have a jump, etc. If we write naively the usual expression for 
the variation of the action in the case of delta-function fluid distribution, then many 
terms in it look meaningless within the theory of distributions (^-functions multiplied 
by discontinuous functions, etc.). Some ingenious calculation of all variations might 
still lead to meaningful expressions. Instead, we resort to a simple trick by which 
the problem is avoided: we fix the spacetime coordinates of the discontinuity surface 
S. In this way, the surface of the discontinuity is formally made to an 'externally 
given' boundary. The fields z"-{x) and g^v{x) will satisfy simple boundary conditions 
at S, and these conditions will be 'inherited' by their variations. Such a strategy is 
possible within the general relativity, because it can be considered as a partial fixing 
of gauge. Indeed, any change of the coordinates of the discontinuity surface S can be 
considered as a superposition of a transformation of coordinates in a neighbourhood 
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of S keeping the physical fields fixed, and a change of the physical fields keeping 
the coordinates fixed; the first step is just a change of gauge. The dynamics of the 
surface is determined by the form of the metric near and at the surface. 

To be more specific about the boundary conditions, let us choose the coordinates 
z"" in Z such that is given hy = 0, and the coordinates in M such that S 
is defined by = 0. Thus, for the matter fields, we require 

z% = 0, 6z^\^ = 0. (49) 

It follows that 

= 0, = (50) 

in the Euler picture. We further assume that 

Condition 1 the spacetime (M, g) is asymptotically fiat and globally hyperbolic; 

Condition 2 the metric gfj,u{x) is piecewise in M, the only discontinuity being 
that its second derivatives jump at S. 

Then the variation 5g^v{x) satisfies analogous Condition 2. 

The total action for our fiuid-gravity system is / = + /g- Here 

lovrG Jv ^ 

where G is the Newton constant and R is the curvature scalar of g^i,. The function 
R{x) can have a step discontinuity at S. is obtained from of Eq. (|]) after 
the following substitution 

z^ix", x^ = 0)= 0, i^(x", = 0) = 0, z\{x'^, x^ = 0) = 0, (51) 

a = 0, 1, 2 and A = 1,2. Thus, 7^ contains less variables than J^. 

The integration volume V is chosen to be bounded by two Cauchy surfaces, 5*1 and 
5*2, and by a timelike surface S"*" (which will be eventually pushed to the infinity). 
Let the coordinates x'^ be adapted also to S"*" so that S"*" is defined by x^ = r~^. The 
matter boundary S divides V into V~ and V~^, and Si into S~ and 5*^^. 



2.2.2 The variational formula 



The variation of the gravity action Ig can be obtained from the following fundamental 
lemma that has been shown in iT^ . 
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Lemma 1 Let the integration volume V of the action Ig be bounded by two spacelike 
surfaces Si and S2, and by a smooth timelike surface S; let be some coordinates 
in V , in Si, in S and rj^ in dSi = S fl Sj. Then 

51, = --^ I d^xG^^'dg^, 



+ 



16n GJv 

167iGJs2 IottGJsi 

If If 
— — / d^ri y/X5a - — — / d^r] y^5a 

/ rf^e7«/35Q"^ (52) 



where 



G>^'':=^\{R>^''-]^g^^R), (53) 

is the Ricci tensor of the metric g^^, qki is the induced metric on Si written 
with respect to the coordinates y^ , q its determinant, 

dx^ dx" 
K = q'^Kki, 

is the future directed unit normal to Si so that K^i is the second fundamental 
form of the surface Si, 'ja/s is the metric induced on S written with respect to the 
coordinates 7 its determinant 

Q"^ := Vl^(^7"^-^"^), (54) 
dx^ dx^ 

Laf3 = riif,.^-^-^, 



/3) 



fh^ is the external (with respect to the volume V) unit normal to S so that Lap is 
the second fundamental form ofE, A is the determinant of the 2-metric Xab induced 
on dSi written with respect to the coordinates rj^ and a is defined by 



a 



-a.Tcsinh{g ^^n^m'^) . 



The Lemma 1 is completely general, independent of the form and description of 
the matter; it determines the 'gravitational part' of the variation formula that we 
are going to derive. For the 'matter part', we can use the formula (0) in which the 
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surface integral along S is left out. Indeed, F{z) = for our special coordinates 
and the boundary condition gives FaSz"" = 0. Hence, 



Jv- V 9z'' ^ dz'^ J Js^ dzS isf dz^ 

+ I f J^xT^'dg,,. (55) 

Our next task is to rewrite the surface integrals in Eq. (^5]) in a covariant way. 
We define, in analogy with Eq. (pTD, 

Pa : 



dz"" 



By a similar argument as in Sec. p.l.3| , pa are surface densities, and the covariant 
forn of the integrals is: 

/ dS^5z''-[ dS^Sz''= [ d'ypjz''- [ (fyPaSzr (56) 
Js- ozq Js- azl Js:^ Js- 

The relations between the old and new matter momenta will play some role. They 
can be summarized as follows: in V~ , we simply have 

Pa = Pa, (57) 

whereas at S, 

PA = PA, P3 = p{PA, 4), (58) 

where p is some function of the variables indicated. Eq. (|57|) and the first Eq. of 
(0) follow directly from the definitions, if the substitutions ( pT] ) is made in the 
expressions on the R.H.S.'s. As Lm does not depend on at the boundary, there is 
wo psis; Psis as given by Eq. (^) with the substitutions (|5T| ) is, however, non-zero 
and it can be expressed as in the second Eq. of (|55|). Let us give a proof. The 
solution of Eq. (^) with respect to reads 

i» = i'^(p.,z^4). 

At S, we must have z^ = 0, so we obtain one constraint for the functions Pals^ 

i3(p,,^«,4) = 0. (59) 



This can be solved for psis; the second Eq. of (|58D is the solution. 

As an example, we work out the explicit form of Eq. (|59D for the dust. We easily 
obtain from Eq. (|l]) 

f = -hX-^x'lz^ 
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if we observe that the intermediately resulting terms can be expressed by means of 
^^-derivatives of the determinant of the matrix z'^. Eq. (^) yields 



h 



n 



X\\9 



-goo + 2NaZ'' - qabz^z^ 



where 



Na 



Nkxl 

k I 



Then, for p = jjn, we have from (PT| ) and 



Pa = -/i 



SO that 



X^\/\g\n 

X\/\9 



-npb 



and the desired constraint (59) reads 



X'^N./q , 



= 0, 



where 
and 



~ab ._ kl a b 
1 -—I Zj^Zi 



N" := N'^zl 
For n, we have to insert from Eq. (|30D : 

1 _X^^ 



ji^h^X 2 + q"-^PaPb- 



(60) 



Thus, the constraint ( [60| ) can be written as a quadratic equation for p3, whose 
general solution is 



P3± = -^Pa±N" 



g33 [a^2^33 _ (]V3) 



where 
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and 

^AB = Qab 

is the metric induced on dS. Only the lower sign is admissible, as we can easily see 
from Eq. (0). Hence, finally. 



J g33 (^]V2g33 _ (iV3)2 



Eqs. (|52D, ( 155D and (|56D imply the final variational formula for our gravity-fiuid 
system: 



51 



+ 



4 / 
a X — — 

V- \ oz'^ 



1 



IQ-kGJv 



ST 



d^yqkiSrr''^ + 



d^yqkiSn 



kl 



IGttG Js2 167rG Jsi 

^ / d'^lapSQ-^ + l f d'xT^'^Sg^ 



167 



(61) 



From the formula (^) we can read off the field equations; within V , we have: 



dz" ~ ^ dzf, 



0, 

SttGT'^^ 



and within , we have: 



G'^^ = 0. 



(62) 
(63) 

(64) 



Apparently, the surface equation (|T8|) has been lost. However, using the boundary 
Condition 2, we easily find that is continuous at S. Hence, 



lim G1 







and the surface field equation follows from Eq. (^). We also observe that the 
dynamics can be completely shifted to the gravity if the ideal fiuid is described by 
co-moving coordinates everywhere in V~ . 

Putting everything together, we obtain in an analogous way as in Sec. [^.1.3| : the 
field equations are equivalent to the following relation (which is an analogon of Eq. 
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(5.16) of H) 

5C = [ d'yipjz^y--^ [ d'yiqki67T''y + -^ [ d'v {VXSa)- 
Js- IdttG Js ottU Jds 

^ ^ (fvio.pSQ''^ (65) 



167rG Jds 

where C is the Lagrangian of the system 



C := / dSL. 

Js 

Eq. ( p^D is the variation formula for our system. 
2.2.3 The Hamiltonian formalism 

Eq. ( |65D is a good starting point for the Legendre transformation to a Hamihonian 
formahsm. We define one of the conceivable total Hamiltonians for our system by 

H:= f d^ypaz^ - / d'yqkiTr'' + ^ / d'v VXa - C. (66) 

Js- 167rG Js SttG Jas 

Then, the field equations can be obtained from the variation formula 

6n = f d'y{z''6p,-pjz^) + -^ f d'y{qMSn''-7T'''6qM) 
Js- 167rG JS 

+ ^ / d'v:;^{cySX-X6a) + -^ [ d'v^af^SQ"^. (67) 
SttG Jas 2v A 167rG Jas 



To find the explicit form of the Hamiltonian we use the lemma 



Lemma 2 In the adapted coordinates defined in Sec. \2.2.% the following identity 
holds at any S = St in the volume V 



I dSqkiT^^^ -2 I ddSy/Xa = 
Js Jas 

- 2 f dS J\g\R'^o - 2 / ddS V^L°, (6J 
Js Jas 



where L"^ is the second fundamental form of the boundary dV corresponding to the 
normal oriented outwards from V , dS = dx^dx'^dx'^ and ddS = dx^dx"^ . 



The derivation of this identity is given in |jT^ (Eq. (6.3)); the form (|68|) is easily 
obtained if one uses the equation Q^^Qab — Q^^Qoo = "which follows from 

Eq. (1541). We also observe that 



/ dS PaZ Cm — "^r. 

Js- 
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where Tim is obtained from Tim as given in the Sec. |2.1.4| by the substitutions 
(0), ( |57D and (|58|). Thus, the substitution for is discontinuous near S. It 
could, therefore, seem that the corresponding Hamiltonian density Hm would not 
be continuous at S, but this is not true. The reason is that -ffmis does not depend 
on p3, if the conditions (^) are satisfied: 



dpa 

Hence, we have from Eq. 



i'ls = 0. 



Ti-m = — I dS Tq. 



Collecting all results, we obtain finally: 

« = //^ H + s^G^S) + ^ das VHL". (69) 

This is the full 'off-shell' Hamiltonian of our system. If Einstein's equations hold, 
its value is just the surface integral. 

Let us rewrite the volume integral in a covariant form. Any tensor density VF^ 
satisfies the identity 

The unit future-oriented normal covector ri^ to S has the components = — A^5° 
with respect to the adapted coordinates; it follows that 

= -^{Nw^^ + N'^wi), (70) 

where 

w^^ := -^W^n^n,, 



1 

1^1 



Thus, the volume integral can be written as 

1 „n\ 1 



where 

C ■= ^(G^'^ - ^TiGT^"')n^n^ (72) 
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and 

Cfc := -^{Gt - 8nGTi:)n^ (73) 
V 1^1 

are the super- Hamiltonian and the supermomentum of our system (or scalar and 
vector constraint functions). 

Let us return to the formula (^), which not only implies the field equations, if 
we perform the variation on the L.H.S. and compare the result with the R.H.S., 
but it also determines the so-called control mode (see, eg. |]13|) and the type of 
boundary value problem for the field equations. We observe that this mode is 
a kind of "curvature-control-mode"; it amounts to keeping fixed (controling) the 
external curvature Q""^ at the boundary 'L'^ (see [|14|). Such a boundary problem 
for Einstein's equations has not been studied. To pass to a more natural and in fact 
conventional approach, we have to perform an additional Legendre transformation 
at the boundary [Q: 



IottG Jds 



so that we have, using also Eq. ([71|), finally 

1 /• ,w„„ . 1 



and 



^ = FTT / + - TFTT / ddS g°°7oo (74) 

SttG Js IottG Jds 

Js- IottG Js 

+ I ddS ia5X - X5a) 

8ttG Jds 

+ TTTTT / rf55(7oo5Q°° + 27oA5Q°^-Q^''57AB). (75) 
IottG Jds 

The last surface integral in Eq. (^) will result in the A.D.M. energy, if the limit 
S"*" — »• cxD is carefully performed (this has been shown in [|14|)- The last one in Eq. 
([75|) defines the way of controle: (5°°, Q^^ and 7ab are kept fixed at the boundary. 

The transformation to the Euler picture in the matter part of the Hamiltonian 
is straightforward; let us denote the resulting Hamiltonian density by H'^. Most 
formulas of Sec. |2.1.5| will result in the analogous formulas for H'^, if the substitutions 
(^), ( 1^ and ( 1^ are performed in them. For example, we have to use the modified 
formula (^), which will read at S: 



Tl = -pBzi, n' = -Pbz^ - p{p, z)zl, 

etc. Only those formulas that contain derivatives of if^ with respect to the variables 
which are not contained in H'^ (like psis and x^ls) need some care to be properly 
transformed. 
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We observe finally that the field equations derived from Eq. ( |75D (or an analogous 
equation of the Euler picture) will have the form of canonical equations. This is 
interesting, because Eq. (^) does not seem to have such form. The variations in 
Eq. 



[5]) (or that of the Euler picture) must satisfy the boundary conditions (^91) 
(or (|50|)) and will, therefore, lead to trivial surface equations; the volume equations 
alone have the canonival form. For example, the would be counterpart of Eq. (^5]) 
originates from the term 



in the variation of H'^. As 6x^\ 
reads 



Further, a counterpart of Eq. 
is equivalent to 



= 0, the only equation which is implied thereby 
dH' 



dx4 



0. 



is valid for this derivative, so the above equation 



0. 



Moreover, z^\-£ = 0, so the equation reduces to T^ls = 0. However, these compo- 
nents of Tj^ vanish identically at S, besause = (cf. Eq. (||)). Thus, there is 
no surface field equation. 



3 Fluid shell 

In this section, we are going to describe the dynamics of a delta-function distribution 
of fluid. The matter will be coupled to the dynamical gravity from the start. We 
shall consider a special case: just one shell in vacuum; a generalization to more shells 
surrounded by a piecewise smooth matter is straightforward as far as the shells do 
not intersect. 

3.1 Action in the volume form 

The shell can be represented as a delta-function singularity in the mole density 
h. The action can then be written as a volume integral of the same form as for 
a regular distribution of matter. This holds also for the gravitational part. We 
shall give a more detailed description of this volume form and then transform it 
to a combinations of volume and surface integrals, where no delta- functions will 
feature. This may be useful, because much more general choice of coordinates is 
then allowed. Indeed, the ^-function method works only if the coordinates are such 
that the corresponding components of the four-metric are continuous. Further, the 
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coordinate position of the shell — the three-surface in the spacetime M and the two- 
surface in the matter space Z — will be kept fixed. Here, everything can be repeated 
what has already been said in Sec. |2.2.1| about this point. 



The matter space Z remains, therefore, three-dimensional first. Let be a two- 
dimensional surface in Z on which the matter is concentrated. Let be coordinates 
adapted to so that the equation = Zq determines S^. Such coordinates are 
determined up to a transformation 

z'^ = z'^{z''), z'q^ = z'\z^,zl), 
where Zq is a constant independent of z^. Then, 

d{z^,z'^,z^) _ d{z^,z^) dz^ 
d{z'\z'^,z'^) ^ d{z'\z'^) ' d^' 

We decompose the molar density h in the adapted coordinates as follows 

h = hs6iz^-zl), (76) 

where hs{z^,z'^) is a two-dimesional density on S^. If we change the adapted coor- 
dinates, we have 

~ d{z\z^,z') d{z\z^) ,^)^'' -I .5f53_~3^ 

SO the decomposition is independent of the choice of adapted coordinates, and de- 
fines, in fact, a two-dimensional matter space with a two-dimesional mole density 
hg] later, we will pass to this space. 

In the spacetime M with coordinates x^, the matter fields are z"'{x'^); the shell 
occupies a three-dimensional surface S, which can be described by the embedding 
functions x'^ = a = 0, 1, 2, or by means of the equation z^{x'^) = Zq. Later, 

we will pass to the the matter fields z^ = z^i^^'^), where z^^^"') = z^{x'^{C,°')). 

The gravitational field is described by the metric g^jj{x)\ we require: 

Condition 1' the spacetime {M,g) is asymptotically flat and globally hyperbolic; 

Condition 2' there are coordinates in a neighbourhood of each point of S such 
that the metric gfj,u{x) is C° everywhere, piecewise C°°, so that the only dis- 
continuity is a jump in the first derivatives at S. 

The second derivatives of the metric will then have a delta-function singularity at 
S so that the Einstein's equations can be satisfied. 
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The total action for the system consisting of the shell and the gravitational field 
can then be written in the following form 

I = -J— f d'^xJ\g\R- I (fx J\g\ne{n). (77) 

One can use this volume form of the action to derive the equations of motion. 
However, there is also a 'surface form' of the action; we are going to derive this one 
in the next section. 



3.2 The surface form of the action 

The coordinates satisfying Condition 2' are not uniquely determined. We can use 
this freedom for the derivation of the surface form; the tool will be the adapted 
spacetime coordinates x^, defined by the property 

x^ = z^x), C = x'^Ie- 
Then, the induced metric 7^/3 on the shell is 

7a/3 = fl'a/3|s (78) 

and its determinant 7 is related to the determinant g of the four-metric g^^, by 

1 = 9- (79) 

3.2.1 Matter action 

Formula (|I]) together with Eq. (^) give 

Hence, = 0, and as 6{z^ — z^) = z^6{x^ — Xq), j" can be written as 

r=mx'-xi). (80) 

The mole density n can be calculated from Eq. (^. We obtain 



^-g^.vj^j" = \/-9a/33fjs6{x^ - Xl) = \/-lat33?3sS{x^ - xl) , 

and Eq. (|79D yields: 

^^/-ic^JsJss/g^Six'' - xl). 



n 

'|7| 
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We define the surface mole density Ug by 

1 

so tliat 

n = ns/g^5{x'^ - xl) , (82) 

wfiere \/g^6{x^ — Xq) is already a scalar witli respect tlie reparametrizations of x^. 
For the velocity u^, we have the expansion 



where 



u^' = w"e^, (83) 



e 



lapV^V^ = -1 



SO that 
and 

Then the transcription of the matter Lagrangian density is straightforward: 



Lm = -^/^ne{n) = -\/\'y\nses{ns)6{x-^ - Xq), 

where es{ns) is the energy per mole of the shell matter. We define the surface 
Lagrange density and the surface mass density by 

Ljn := Ls6{x^ - Xq), Ps := nsCsins), 

so that 



Ls = -^lilpsiris). (84) 

If we perform the trivial integration over x^ in Im, the matter action becomes a 
surface integral 

Im = - d^^ \/W\r^ses{ns)- (85) 

This expression is invariant under the transformation of coordinates at the shell. 
The action Im can be varied with respect to the shell metric jafi with the result 

On \ 
llp'sjr-^ hap- 



In analogy with the formula (|^), we obtain 

dus 1 



d^ap 2 
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-i.nsiio.p + v'^v^), (86) 



and we have 

5Im = -jj^iTf5^^p. (87) 
where is the surface stress-energy tensor, 

Tf ■■= [psv'^v^ - a{^o.p + v'^v^)) , (88) 

and 

■= nsp's - Ps = -nle'^ 
is the surface tension (negative two-dimensional presure). Moreover, it holds that 

We obtain easily relations analogous to the Eqs. (|11|), ( [T^ ) and (|13|): 

r-(0 ^ 2^ (89) 



(this is a form of Eq. (pS])); 



g^--,^ + 2^-„, = W?, (90) 



because Lg is a three-density on S, and the Noether identity 

dL 



= Ls5t - ^^zl (91) 



3.2.2 Gravitation action 



The next task is to rewrite the shell part of the gravitational action Ig in the surface 
form. The following lemma is vital. 

Lemma 3 In the adapted coordinates that satisfy Condition 2', the delta-function 
part of the gravitational Lagrange density is given by 

g\R = --^J\^\[L]Six' -xl) + ..., (92) 



where the dots represent regular terms, L = '-y^pL"-^ , L""^ is the second fundamental 
form ofT, corresponding to the normal oriented outwards ofV~ and the abbreviation 
[/] := /+ — /_ for the jump of a quantity f accross the shell is used. 
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The proof of the lemma is relegated to the Appendix. Eq. ( ]92| ) implies immediately 
that the gravitational action can be transformed to 



h = I ^ ^\R - ^ / rf'e V^W- (93) 

IoTTG Jv+yjV- QTlKj JS * 

Eqs. (^5]) and ( ]93|) give the total action in the surface form: 



- ITF / \/H W - / ''"^ \/b\n,e,M. (94) 

This action functional is equivalent to that given by Eq. ([77D, if the coordinates 
satisfy Condition 2'. It has, however, two advantages in comparison with ([T^): 

1. all integrands in are smooth, 

2. it is valid and can be used with more general coordinates, namely arbitrary 
smooth coordinates x± within and arbitrary coordinates within S. 

The fields in the action (0) are the matter fields z^(y) on S (observe that the 
fictitious field disappeared from the action), the gravity fields g^u{x) in and 
la(i{y) in S. The metric has to satisfy the so-called continuity relations 

7«/3(0=(v|^|^(-(0)j = [9%-^-^i-.{0)) , (95) 

where the symbols ()^ denote the limits from the volumes towards S. The role 
of the continuity relations ( P^D is to define the configuration space of our system 
similarly as a control mode or some fall-off conditions do. The embedding functions 
^±(0 cire fixed; their variation is zero. 

We also have to specify the integration volumes; this will be done in analogy to 
Sec. |2.2.1| : the volume V is chosen to be bounded by two Cauchy surfaces, 5*1 and 5*2, 



and by a timelike surface E"*" (which will be eventually pushed to the infinity); the 
surface S separates V in two parts, V^, and the surfaces Si into Sf] the intersections 
of S with Si will be denoted by 9Sj and we will assume that they together form the 
complete boundary of S; the intersections of S+ with Si will be denoted by dSi and 
we will assume that they form the complete boundary of S+. 

This form of the action will be our starting point to the derivation of the field 
equations as well as the Hamiltonian formalism. 
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3.3 The variational formula 

The variation of the matter part Im of the action (p4|), if we calculate in the coor- 
dinates that are adapted to the surfaces dT,i by = U and rj"^ = at SSj, 
is 



JdT.2 OZq JdT,i OZq 



dLs r .0 dL 

We define the matter momenta pa by 



dL. 

As Pa is a well-defined two-surface density (cf. the discussion below Eq. (PT|)), we 
obtain the covariant formula: 

+ / d^7]pASz^ - I d^r]pASz^. (96) 

JdT.2 JdT,i 



To calculate the variation of the gravitational part Ig of the action (p^), we first 
rewrite the surface integral in Ig with the help of the trace part of the Eq. 



Q = yh\L 

as 



Xd'^v1^[Ll = -^Xd'«,„,|Q°a (97) 



SnGJs 
so that 

Then, we apply Lemma 1 (Eqs. ( ^2] ) and (|5^)) and Eq. ( pS]) with the result: 
lovrG jy~ lovrG jy+ 



IGttG Js 

IbTrG J52 167rG Jsi 

f d^T] X6[a] + f d^r]\6[a] 



+ 

»7r 



^ / ci^r^ A5a - / rf'r/ A5a - / lapSQ^^ ■ (99) 

tG ovrG JdSi IdttG js 
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Here, qki{y) and Squiiy) are continuous along Si, but n'^^y) and 6Tr''\y) have a jump 
at (9Sj. Eqs. ( ]96| ) and ( pQ] ) imply the following formula for the total action: 



" = -liSG t''*^ -lebX.^'^ 



SttG 7as2 SttG Jasi 

+ ^ f d'v ^Sa - ^ / d'v V\5a - / d'i la^SQ^^ ■ (100) 

SttG Ja52 ovrG Jasi IdttG ./s 

From the variational formula (|100| ), we can read off the field equations: in V~ 
and V+, 

G^"" = 0; (101) 

at E, we obtain 

[g"/3] = 87rGTf'5. (103) 

Eq. ( |103| ) is the well-known dynamical equation for thin shells |jlT|; we shall refer 
to it as 'Israel's equation'. It may be considered as the singular part of Einstein 
equations, corresponding to S-\ike sources. Eq. (|102|) can be interpreted as a three- 
dimensional stress-energy conservation: we can show in an analogous way as in Sec. 



2X^ that 

where Vq, is a covariant derivative associated with the metric 7q,^. This identity 
implies that Eq. ( |102[ ) is equivalent to 

V.T° = 

(which comprises only two independent equations). 

Eq. ( |100| ) implies a generating formula for the field equations analogous to Eq. 
(p^). In order to derive this formula, we first have to introduce a foliation of the 
integration volume in (|100| ). This is an arbitrary smooth family of spacelike surfaces 
St such that Si = St^ ; we allow for the surfaces St having a cusp at dJ^t = S fl S't so 
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that the normal can have a step discontinuity there. This leads to jumps in vr'^' 
and a across dT,f 

We also have to introduce adapted cordinates x'^ so that the surfaces St are given 
by a;° = t, S by a;^ = and by = r+; further, = x^\st, = a^"|s, 

= = x^\dT,t ^iid 7]^ = x^ldSf Observe that the full four- metric need 

not be continuous across S with respect to these coordinates. Observe that t = 
const is a continuous surface intersecting S and S+, and d/dt is a continuous vector 
field everywhere. 

Eq. (H) and 

imply for C: 



rt2 

/ dtC 
Jti 



167rG ^ IbvrG Js+ ^ 

^— I d^r]J\^\[L]- f d^r]J\^\n,esin,). (104) 



»7r 

If we rewrite Eq. ( |100| ) in the form 



61= r dt6C. 



we obtain the variation formula for our system: 

167rG J St SttG JdY., 



IGvrG JSt SttG JdT,t 

Jast SvrG JdSt 



1 



p , rf'r/7a/55Q"^- (105) 
IbvrG Jast 

Performing the variation in (|104|) and comparing the result with the R.H.S. of (|105|) 



recovers the definition of momenta and the field equations. 

3.4 The Legendre transformation 

Let us define the Hamiltonian in analogous way to Sec. 2.2.3| : 



n = -C- / d^yq,i7^'' --^ [ d'r^ v^[«] 

167rG Js SvrG Jsns 

+ d^r]y/\a+ f d^r]pAZ^, (106) 

where [a] := — a+ — a_ at S, a;+ is defined by the normal to S that is outward to 
and the future normal to S~ , is defined by the normal to S that is outward 
to and the future normal to S~^. 



33 



To calculate the variation of Id, we have to regroup terms in Eq. ( |105| ) 
1 



5C 



IGttG 
1 

IGttG 
+ 2 



+ 



Js- Js-nT. 
Js+ JS+ns 



+ 



A\. 



5ns 



1 



167rG Jds 



I d^'qj^fsSQ'^^ 
JdS 



Then we vary Eq. ( |106D , substitute for 6C the regrouped expression, and apply 
Lemma 1 to each of the two volumes V^; the result is 



sn 



1 



levrG Js 



kl\ 



+ 



JsnT. IdttG JsnE+ A 



a- 



6X, 



(107) 



IGttG J5ns+ 

This is the Hamiltonian variation formula for the field equations. 

Let us compute the value of the Hamiltonian. To this aim, we insert from Eqs. 
(0) and (|97|) into Eq. ( |106|) and apply Lemma 2. A simple calculation leads to 



H 



SttG Js 



d'yG',+ 



1 



SttG Js 



d'yG', 



1 



SttG JsnT, 



JsnT, SttG Jsr\T,+ 



(108) 



The term ~^^[Qo] interesting interpretation. It has been obtained as 

the sum 



1 



1 



87rG^^°^ 167rG^~"^ ' SttG 



-IQ] 



1 



The first summand is the delta-function term in the gravitational Lagrangian density 
in the volume form, 

IottG ^ 

and the meaning of the second is given by the following lemma: 

Lemma 4 // the four-metric is continuous in the adapted coordinates, then we have 
atT,: 



8nG 



\9\Ro 



1 



SttG 



[Q S{x' - xl) + 



(109) 



where the dots represent regular terms. 

Thus, the term is the delta-function part of the expression 

SttG^o 
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(cf. (p3D). If the assumptions of Lemma 4 are satisfied, then the three first integrals 
on the R.H.S. of ( |107| ) can be written in volume form just as 



1 



SnGJs 



The Legendre transformation at S+ similar to that in Sec. |2.2.3| can be performed 



exactly as in Sec. p.2.3| . The transformation to the Euler picture in the matter part 
of the shell Hamiltonian ( |106| ) is much simpler than the analogous transformation 
of the step Hamiltonian (|66|), because all formulas of Sec. |2.1.5| remain valid, they 
must only be rewritten in three spacetime and two matter space dimensions. 



3.5 The form of the Hamiltonian 

In this section, the Hamiltonian ( p.08| ) will be expressed as a functional of the canon- 
ical variables A^, A^^, qki, t^^\ A and a. 

Observe that the formulas (107) and (|108|) are valid in any coordinates that are 



adapted to the foliation and that make the embedding formulas for E and S"*" time 
independent. More specifically, the coordinate t must be constant along the surfaces 
S, the embedding formulas for E and S"*" must read 

t = e. yi = yi{e). 

and the embedding formulas for S* fl S in S is 

e = const, e = V""- 

Let us recall that one important point of our method is that the boundaries are time 
independent in the above sense and their variations are zero. 

The 2+1 decomposition of the metric at S and S"^ is analogous to that of 
Qfj^iy. In particular, we define the (surface) lapse u and the (surface) shift uk by 

/ 1 







,2 1 



SO that 

7oo = -i^"^ + X^^i^ki^l; 
the 2+1 deconposition of the continuity relations (pS]) reads 



1. = ^Ni-{N^y, (110) 
UK = N^ei^, (111) 
Xkl = QueiKeiL, (112) 
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where 



and m\. is the unit normal vector to S fl S* tangent to S and oriented from S~ to 
5*+. From this definition, it follows that 

Using the decomposition (|7DD, we can write for the integrand of the volume terms 
in Eq. ([TOSl) : 



9\ ^\9\ 

observe that the R. H. S. is invariant with respect to transformations of coordinates 
y^. The form of G-^-^ and is well-known (cf. 0) 



— 4Vg 2 ' 



'-^k — ~^k\ii 
\9\ 

where R^^^ is the curvature scalar of the metric qki- 

Within S and S+, an analogous decomposition yields for the surface terms: 

n\ 



and 

[Qo] = -^HQ'-^] + ^''{Qk])- 

n\ 



r°Q can be expressed by means of the canonical variables z^, pa, ^kl, ^ and uk 
in a way parallel to Sec. |2.1.4| : the formulas are independent of the dimension of 



spacetime. In general, the form is only implicit. However, the dependence of T°q on 
\kl, ^ and vk can be inferred from the relation 

Tf = -2^t3l (113) 



analogous to Eq. (|27|) . It follows that 

=T,^^, (114) 



dv 
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and 



^^-^T.i. (U5) 



In particular, -^^T^^ and -^Tji- are both independent of v and vk- 



In an analogous way as in Sec. 2.1.4 we obtain easily 



and 

For a dust shell, we have 



1 _^ , X'^^Z^zfpAPB 



and 



Eq. (|116|) then yields 



PAPBZ^zf 

J-sKL 



^fi'if.y + x^^z^zEpAPB 

The form of [Qq] can be given explicitly. Let us observe (cf. [T^) that the normals 
m, n and m are related by 

■fh^ = sinh a + cosh a, 
n'^ = n'^ cosh a + sinh a; 

recall that n is the normal to S in M, m is the normal to Sfl S in (m is orthogonal 
to n), n is the normal to SflS in S and fh is the normal to E in M (n is orthogonal to 
fh). A simple calculation then confirms that the corresponding second fundamental 
forms Kki and Ikl of S in M, S* in M and S fl S' in 5, respectively, satisfy the 
relations 

La/se^e^L = -Kkie^e^L sinh a + Ikl cosh a, 

It follows that ^ ^ 

Q"*""*" = — TT~^~^ sinh a — / cosh a. 
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Qk l=^K = ^,K- 



_ 71 ^ 
Hence, 

go = -z/(\/A7r^^sinha - \f\l cosha) - u^i^/Xn^ + VAa,^), (117) 

where 

depends only on q^i and its first derivatives. We use also the abbreviations 

-yj-fe^ ^fei TT^^ 

where 



Finally, the complete Hamiltonian reads 

1 I" 

— ——^ / (i^?7\/A (i/[7f"'""'" sinha — Z coshct] + i/^[7f^ + a ) 

- / d'7^T^,+ -^[ (fr]Ll. (118) 

The surface term at E"*" is left unchanged; it has to be transformed according to the 
control mode used and/or shifted to infinity. 

The surface super-Hamiltonian T-Lg and the surface supermomentum T-Lsk at the 
shell are given by 

— — „ [tt"'""'" sinh a — I cosh a] + T^"*" , 



here 



The = -^TffiaepK. 
7 
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The geometric meaning of the gravitational part of Tig and TisK can be inferred 
from Eq. (M): 



Q 



±± 



Qk t ~a 



In particular, Lape^e^ is the second fundamental form of the two-surface S'flE corre- 
sponding to the normal m (each two-surface has two independent second fundamen- 
tal forms in the spacetime); hence, the gravity part of the surface super-Hamiltonian 
is the jump in the (two-)trace of this form. 



3.6 Equations of motion 



In this subsection, we calculate the variation of the Hamiltonian ( |118| ) explicitly. In 



this way, we can check if our method leads to the well-known equations of motion; 
moreover, we can study the structure of the canonical equations and constraints at 
the surface S. In varying the Hamiltonian, we must carefully deal with boundary 
terms. 

The variation of the volume integrands can be given the form: 



5(26-°) = CkSN'' + C5N + a^'dq^i + fefc^^vr^' + ^B^^ 



where 



2 



^kl (c\„k „lm „„kl „mn„ kl 



kl 



and 



bki = — (27rfei - nqki) - Cj^qki, 



= — —(N''n^'' + N^7r^''-N'^7r''^)6qki-—Ni67i^'^ 

- N{q''6ri, - q'^Sri,) + ^Nsiq^'q^^ + q^'q^' - 2q^^q'')6qki. 
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Here, jCj^ is the Lie derivative with respect to the vector field X. The comparison 
with the volume term of Eq. ( |107| ) yields the well-known canonical form of Einstein 
equations (cf. ||^): 

Cfc = 0, C = 0, 
tt'^' = -a''\ qui = hi. (119) 
The divergence term contributes to the variation of the surface term at S" fl S by 



-1-/ d'vVXi-B 



t + B'lml 



The following identity can be easily derived 

q''mr{STli - 61TI) = 261 + 21''H\kl - A^^(m^e^5gH)||i^, (120) 

where the symbol "||" denotes the covariant derivative associated with the metric 
\kl on 5*0 2. Using Eq. ( |120|) and the continuity relations (|111|) and ( |112|) , we can 
rewrite the surface term in the 2+1 form 

VXB'^mk = -VX{2n^^N^ + n^^UK){m''m^6qki) - 2^^^^u^ {e)^m^5qki) 
— vA(7r z/+7r u — ly it + tt iVA +n z/jv/A 

+ 2^/XN5l-2^/X{N^mkmi + VKe^mi)5^^K (121) 
The variation of the surface term [Qq] can be written in the following way 

5Ql = 5v I -4^g^^ I - (5z/^ i^Q]A + Rest. (122) 




171 

Similarly, using Eqs. (|114|) , (|115|) and (|116|) , we obtain that 



71 / VVl7 

Comparing the first two terms in Eqs. ( p.21| ) and ( |123| ), we obtain the first three 
Israel's equations: 

[Q^^] = SvrGT,^^, [Qi] = SvrGTi,. (124) 

The last two terms in Eq. (|123|) , if compared with the corresponding surface term 
in Eq. (|107| ) yield the dynamical equations for matter inside the three-dimensional 
spacetime of the shell surface S: 
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'A 



(126) 



What remains from Eq. (|107|) can be written as follows 



[2Rest - VXB' 



rrik 



STrGTf ^ - v^[a]A^^) 6Xkl + VAA^^Ai^L^H. (127) 



A somewhat lenghty calculation starting with Eqs. ( |121| ) and (|117|) gives 

2Rest - ^/XE'^nik = 

-2\/A7f^^(i/sinha - N^)m^m^6qki - \^{n^^X^^{iysinha - N^) 
+n^^N^ + F^N - ulX^^ cosh a + u^^X^^a^M - N ^km^X^^}5XKL 
+2v^(z/cosha - N)5l - 2v^(z/sinha - N^)mkmi6TT''^ 
+2-\/A(— z^vT"'""'" cosh a + ul sinha + z/|j^)5a. 

Substituting this into Eq. ( p.271) , we obtain immediately 

u sinh a± = , u cosh a± = N± , 

[— TT^"*" cosh a + I sinh a] = 0. 
The remaining equations, simplified by (|128|) and (|129|) , read 

and 

X = -2X{Nn^^ - N^I-u^k). (131) 
From the point of view of physical (or geometrical) content, Eq. ( p.28| ) just re- 



KL 



(128) 
(129) 

(130) 



produces the definition of a and it is compatible with the continuity relation (|110|) . 
The three Eqs. ( |13CI| ) are equivalent to the remaining three Israel's equations. Fi- 
nally, Eqs. ( |129| ) and ( |131| ) follow from the continuity relations ( |110| )-( p32D and the 
equation of motion ( |119|) . Indeed, taking hmit of the second Eq. ( |119| ) from both 
sides towards the shell and projecting the result by e^e^, we obtain 



Xkl = ( 2N{nKL - ^^Xkl) + Nk\\l + N^k + 21klN^ 



± 



the expression in the brackets on the R.H.S. must be continuous, hence 

[-N{nXKL - 2nKL) + 21klN^] = 0. 
The trace of Eq. (^) is Eq. QT^ ) and the trace of Eq. ( |T^ ) is Eq. QT^ . 



(132) 



(133) 
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Eqs. ( p.24| ), ( |133| ) and the tracefree part of Eq. ( |130| ) can be written in a more 



symmetric form 

[-KKLsmha + IxLCOsha] = -8nG{fsKL-^T\KL), (134) 

[i^xL cosh a — sinh a] = 0, (135) 

i^K + c^,k] = -SvrGT,-^, (136) 

where K^l = —t^kl + (l/2)vrAxL- They give the jumps of the two independent 
second fundamental forms of the shell 2-surface in the spacetime, one corresponding 
to the normal n in the direction of the shell motion (continuous), the other to m, 
which is perpendicular to the direction of motion. 

From the point of view of the theory of constraint systems, the Eqs. ( [124| ), (|129|) , 
and the tracefree part of Eq. (|13CI|) are constraints. The trace of Eq. ( |13CI| ), 

[a] = AtiGuT^'-Xkl - \ {-N^J^'^'^Xkl + NI + 2Nkm'' + 2z/^a,j^] , (137) 
Eq. (PID and Eqs. (PSQ , (PB| ) are canonical equations. Finally, Eqs. (|TTI1|)-([TT^ 



and ( p.28| ) are defining equations of the Hamiltonian system, analogous to fall-off 
conditions or control conditions. 

It seems that some of the constraints are second class. For example, Eq. ( |129|) 



follows from the variation with respect to a = (1/2) (a+ + a_), which is a Lagrange 
multiplier. Eq. ( |129|) contains this Lagrange multiplier; thus, its Poisson bracket 
with TTa, which is the momentum conjugate to a, and which is also constrained to 
vanish, is not zero (if we extend the system by this momentum). 

Another important observation is that the L. H. S. of Eq. ( |129| ) can be smeared 
only by a function of two variables, because the domain of definition of the L. H. 
S. is the shell surface. On the other hand, Eq. ( |129| ) contains so-called volume 
quantities, namely vr^' and l^^; a derivative with respect to these variables and the 
Poisson brackets of these variables result in three-dimensional (5-functions. Thus, the 
L. H. S. of Eq, (|129|) cannot be smeared so that it becomes a differentiable function 
on the phase space. We call such constraints singular. The best way of tackling this 
constraint may be to solve it for a and insert the solution back into the action (cf. 



T^). A similar procedure exists hopefully for the two constraints which result from 



the tracefree part of Eq. ( |13(J| ) : 

fr^^ - Ivf^^^AM^vA^^) cosh a - (/^^ - 1/^^^Am;vA^^) sinha 

they together with Eq. (|129|) exhaust the singular constraints of our system. It is 
interesting to observe that the total Hamiltonian is a differentiable function. Al- 
though the surface integrals in the Hamiltonian represent two-dimensional smearing 
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of some volume quantities, the presence of volume terms and the continuity relations 
between the surface {u, u^) and volume {N, N'') smearing functions guarantee an 
effectively three-dimensional smearing of all volume quantities. 

These difficulties and the related problem of the Bergmann-Dirac analysis |17| 
seem to be non-trivial; we will try to tackle them in a future paper. 
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A Proof of Lemma 3 

Let us calculate the delta-function terms in the expression Lg = h^^^Rfj,^, where h^'^ 
is defined by 

h^" ■■= TT^I^^rV^ (138) 
lovrG 

For this aim, we have to isolate the second derivative terms. A simple calculation 
gives 

h^^'R.u = d,{h^^Ap + . . . , 

where A'^^ is defined by 

^tu ■= - StK)K- (139) 

and the dots represent regular terms. Thus, the delta-function term in h^^^R^y is 
given by 

W"[AlM^'-xl). 
Thus, we are to prove the following identity 

h''[AU = -^^\[Ll (140) 

where Lap is the second fundamental form of S and L = 'y'^^Lap- Observe that the 
R.H.S. of (|140|) is written in a thre-covariant form. From the definition of ( |139| ) of 
the quantity A^^, it follows that 

^33 = -rL, (141) 
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1 

2 



(142) 
(143) 



Let denote the unit normal vector to S oriented in the direction of inceasing 
that is outward with respect to V~: 



1 



St 



i33 M 



Then, 



m 



SO 



(144) 

(145) 
(146) 



h^^ is a tensor density, hence 

but it is a tensor density formed from the components of the metric tensor, thus 
Vr/i^'^ = for all k, fi, v. Setting /i = 3, z/ = K = ain this equation and using Eqs. 
([141|), ( |145|) and (|14|), we obtain that 



^33 



Similarly 



A 



a3 



2/^33 



(147) 



The metric Qfiuix) is continuous across the shell, so will be h^'^ and the tangential 
derivatives of h^^'^. It follows that 



[^33] 



9 



a(3 



[Lap] , 



[A. 



9 



a3j 



'9 
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[Lap], 



(148) 



(149) 



where we also have substituted for h'^'^ from ( [L38| ). Moreover, using Eq. (|146|) , we 
have that 



[Al 



9 
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The Eqs. (|4|), (1T4|) and (|150]) imply 

1 



(150) 



h^^[AlJ\ 



SvrG 



3a 3/3^ 
33 I _ ^ ^ 

33 



^ir" I 9 
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However, the following well-known relation holds: 

^33 

from which, if Eq. ([79|) is used, the identity (|140|) follows immediately. 



B Proof of Lemma 4 

If we rewrite Rq in terms of the connection, 

Ro = a'' {d,T;, - dov;^ + ... 

where the dots represent terms that do not contain second derivatives of the metric, 
we obtain immediately that 

R'o = 9'Tlo]Six'-xl) + .... (151) 

We easily find: 

The following two equations are easily verified: 

La/3 = (^-^^lalB — Qafi^ ■ 

Then, some computation leads to 



7 



Finally, using Eqs. ( |75|) and (131), we obtain Eq. ( p,09| ) immediately. 
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